Mesoscopic phases of dipolar ensembles with polar molecules and Rydberg atoms 
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We discuss the realization of mesoscopic phases of dipolar gases relevant to current experiments 
with cold polar molecules and Rydberg atoms confined to two dimensions. We predict the existence 
of superfiuid clusters, mesoscopic supersolids, and crystals for a small number of trapped particles, 
with no counterpart in the homogeneous situation. For certain strengths of the dipole-dipole inter- 
actions, the stabilization of purely non-classical crystals by quantum fluctuations is possible. We 
propose a magnification scheme to detect the spatial structure of these crystalline phases. 

PACS numbers: 



One of the most remarkable achievements in recent 
AMO experiments is the preparation of cold ensembles 
of polar molecules in the electronic and rovibrational 
ground state P, The distinctive feature of polar 
molecules is their comparatively large electric dipole mo- 
ments (of up to a few Debye) implying strong dipolar 
interactions. These interactions can be manipulated via 
external DC and AC microwave fields by coupling to ex- 
cited rotational levels of the molecule. In combination 
with reduced trapping geometries this opens the door to 
a study of strongly correlated quantum phases with de- 
signed long range interactions, e.g., for Bose systems the 
superfluid-crystal quantum phase transition in 2D as a 
function of an induced dipole moment While ground 
state Alkali atoms exhibit much weaker magnetic dipole 
moments huge electric dipole moments d ^ in? are 
present in their high- lying Rydberg states, where n is 
the principal quantum number Thus, below we 

propose that strong correlations can be observed in an 
atomic gas of ground state atoms by weakly admixing 
with laser light these Rydberg states. In contrast to po- 
lar molecules, these dressed atomic gases exhibit decoher- 
ence and heating mechanisms from spontaneous emission 
(r ~ and inelastic collisions. However, for strongly 

correlated phases, e.g. in the crystalline phase, the char- 
acter of the many body wavcfunction can strongly sup- 
press the probability for close-cncountcr, possibly harm- 
ful, inelastic collisions, and thus stabilize these phases. 

In light of such new experimental developments, we 
study in this letter mesoscopic quantum phases and 
crossover between these phases in small clusters of dipo- 
lar particles confined to 2D, in the presence of an in-plane 
harmonic confinement [see Fig.[T]Ja)]. By means of Quan- 
tum Monte Carlo simulations based on the continuous- 
space Worm Algorithm we show the existence of 
quantum phases which have no counterpart in the ho- 
mogeneous situation. These include superfiuid clusters, 
mesoscopic supersolids, ring-shaped crystals, and classi- 
cal crystals with unusual rotational symmetries. We find 
that for certain particle numbers and intermediate inter- 



(b) 






T=\ 



;r/r„ I ^ 









•7-=2.5 



■1 'i I x/r„ I 





FIG. 1: (color online) (a) Experimental setup (scheme): po- 
lar molecules or Rydberg-dressed atoms are confined to 2D 
by a strong confining laser beam, with dipoles polarized per- 
pendicular to the plane. In-plane harmonic confinement is 
provided by, e.g., the beam waist, (b-e) Monte Carlo snap- 
shots of the density of particles in all mesoscopic phases for 
= 13 dipoles, as a function of the effective mass r. (b) 
superfiuid; (c) supersolid; (d-e) ring-like crystals; (f) classical 
crystal. 



action strength quantum fluctuation can stabilize purely 
non-classical crystals. Finally, we propose schemes to de- 
tect these crystalline phases. 

We consider a setup where TV dipolar bosonic particles 
are confined to a 2D_plane by applying a strong trans- 
verse trapping field 0, S|, e.g a ID optical lattice, and 
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are aligned perpendicular to the plane, with a DC in- 
duced dipole moment d = ^/D. We assume an additional 
in-plane parabolic trap with frequency co, as realized by 
a magnetic dipole trap, or a single site of a large spac- 
ing optical lattice. For pure dipolar interactions the 2D 
Hamiltonian for N particles with mass m is 
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which is the sum of kinetic energy, parabolic confinement, 
and dipole-dipole interactions. For polar molecules, 
the validity of Hamiltonian ([T]) has been discussed in 
Refs. 0, Q; we will return to a discussion of its valid- 
ity for Rydberg atoms and of relevant experimental pa- 
rameters towards the end of the paper. Defining length 
and energy scales rg = (Z?/tow^)^/^ and eo = mw^rg = 
D/vq = {m^uj^D^y/^ , respectively, Eq. ([T]) can be recast 
in dimcnsionless form as 
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where t = eo/fuj ^ {ro/()^ ~ (mD/fi,^^)^/^ characterizes 
the strength of the dipole-dipole interactions in the trap 
with £ = ^Jh'l muj the harmonic oscillator length. Equa- 
tion ([2]) shows that r plays the role of an effective mass 
representing a control parameter, which can be increased 
by increasing the strength of dipole-dipole interactions, 
or by compressing the trap. 

All the results shown here correspond to a low enough 
temperature T, that they can be regarded essentially 
ground state estimates. Quantitatively, this means work- 
ing at a temperature (much) smaller than the charac- 
teristic energy scale eo. For small r < 1. we expect 
the kinetic energy to dominate, and the cluster to be 
in a weakly interacting superfluid phase. In the limit of 
strong interactions (or equivalently, large particle mass) 
T ^ 1, kinetic energy becomes negligible, and the system 
ground state resembles the classical lowest-energy crys- 
talline configuration obtained by minimizing the last two 
terms of Eq. ([2]). 

We obtain an estimate of the critical for the 
crossover from the superfluid to the crystal by noting 
that for a homogeneous system the superfluid-crystal 
transition occurs at tqm = Dm/h^a = 18 ± 4 [1, Q, 
where rqM represents the ratio of the dipolar interactions 
D/a'^ to the kinetic energy jma with a the mean in- 
terparticle distance. By rewriting Tc = (rgMi/^)^^^, and 
approximating a ^ ^, we obtain the prediction Tc ~ 3, 
which is essentially N- independent. Figure [TJb-f) shows 
our results for r ranging between 1 and 20 and = 13 
particles, which we chose since we found that it displays 
all general features of mesoscopic clusters with N < 30. 
Panels (b-f) are snapshots of the particle density ob- 
tained by tracing over the worldline imaginary-time 




FIG. 2: (color online) Radial density profiles g{r) for the cases 
of Fig.[ljb-f). Inset: statistics of computed particle exchanges 
/ as a function of the number of particles participating to the 
exchange n, for the cases of finite superfluid fraction t — 1 
and 2.5. The supersolid phase with r — 2.5 has finite den- 
sity modulation (figure) and the same /-distribution of the 
superfluid with t — 1 (inset). 



evolution in the Monte-Carlo simulations. Consistently 
with our model estimate, panel (b) shows that for r = 1 
the system is in a superfluid phase, characterized by a 
flat, featureless, density profile, and superfluid fraction 
close to 100%. The corresponding radial density profile 
g{r) is shown in Fig. [H For t > 3 [panels (d-f)] the 
overlap between the various particle probability clouds 
rapidly drops to zero. In particular, for r = 20 we find 
essentially no overlap, as particles are arranged so as to 
mimic the classical crystalline configuration for iV = 13, 
with 4 particles at the center, and 9 outside. In this case. 
Fig. [2] shows two large peaks at the radial position r 
corresponding to the classical equilibrium configuration 
for the crystal, and no density in between the peaks [loj . 

The competition between interactions and confinement 
for intermediate values of r allows for the existence of 
mesoscopic phases with no counterpart in the homoge- 
neous situation. We find that for dipolar particles these 
are: i) a supersolid phase with superfluid fraction one, 
and finite density modulation for t < Tc, and ii) ring- 
shaped crystals with no detectable superfluid fraction for 
T ^ Tc- We have checked that these mesoscopic phases, 
and in particular the supersolid, occur for all particle 
numbers < 30. 

The emergence of supersolid-type behavior for t < Tc 
for A^ = 13 is signaled by a density modulation on top 
of the superfluid background in Fig. [Ijc). This is shown 
quantitatively in Fig. [21 where the radial density proflle is 
clearly more structured compared to the case r = 1. We 
find that the superfiuid properties are unaltered by the 
increased strength of interactions. In fact, the superfluid 
fraction for r = 2.5 approaches unity in the T ^ limit, 
and, as shown in the inset of Fig. ^ no difference is ob- 
served in the statistics of exchange cycles for r 1 and 
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TABLE I: Parameters for possible non-classical/classical crys- 
tal transitions. Columns one to six indicate the number of 
particles A^, the label s for the classical (0) and non-classical 
(1) crystal configuration, the dimensionless quantities /ijv,s 
and !/jv,s of Eq. @, the transition parameter (see text), 
and the ring-configuration, respectively. 



2.5, computed at corresponding low temperatures (i.e., 
same fractions of eq)- 

For T > Tc we observe crystallization of the 
atomic/molecular cloud, consistently with the estimates 
above. This is signaled by a drop of the superfluid frac- 
tion above r = 3. We find that crystallization in these 
finite, inhomogeneous, systems proceeds through inter- 
mediate crystal-like phases, which are characterized by 
the arrangement of particles in concentric rings with a 
fixed number of particles per ring 11| . For small enough 
temperatures, these rings are free to rotate independently 
of each other. Snapshots of these ring-shaped phases are 
shown in Fig.[TJd-e) for t ~ 5 and 10, respectively. These 
two cases differ only for the amount of localization of the 
particle wave- function in the rings. They both display 
no superfluid properties [lo| . 

The crystal-like ring-configurations above character- 
ize the generic crossover between superfluid and classical 
crystal in mesoscopic clusters with < 30. In the follow- 
ing we focus on a few special A'^, for which the crossover is 
characterized by the presence of a mesoscopic analog of a 
first-order transition between a non-classical crystalline 
ground-state and the classical one. These crystals differ 
with respect to rotational symmetry. To our knowledge, 
this is the first demonstration of the existence of these 
phases in mesoscopic clusters. 

In general, for a given A^ there can be several low- 
energy classical equilibrium configurations, which differ 
from the ground-state in the number of particles belong- 
ing to each ring [l3 | . These low-energy configurations are 
found, e.g.; using a nonlinear unconstrained minimization 
of the potential and interaction energy terms in Eq. ^ 
based on random initial configurations Pj. For each con- 
figuration, low-energy excitations can be obtained semi- 
classically by computing the phonons. corresponding to 
small vibrations around the mean particle positions. For 
a particular configuration s (s = and 1 are the classi- 
cal and non-classical configurations, respectively) and A^ 
particles, the mean-energy per particle Ejy^s reads 




FIG. 3: (color online) Non-classical vs classical crystals for 
— 12. (a) Energy as a function of r as computed with 
Eq. Q for the configurations (4,4,4) (red dashed line) and 
(3,3,6) (blue continuous line), corresponding to the non- 
classical and the classical crystals, respectively, (b) and (c): 
Monte Carlo snapshots of the groundstate configuration for 
T = 7.5 < r« and t = 20 > t*, respectively. 



Here, fJ.N,s is the dimensionless classical potential energy 
per particle, while i^n.s — "-^fc/^"^ gives the 

ratio of mean quantum fluctuations to trap energy 
spacing. While the classical ground-state configura- 
tion has always the lowest energy for r 3> 1, it is in 
principle possible that for certain N and intermediate 
T, quantum fluctuations stabilize as a ground-state a 
crystalline configuration (e.g., no superfluidity) with a 
rotational symmetry which is different from the classical 
ground-state. A semiclassical estimate for this transition 
is obtained by the relation = (i^i — h'o)/{^o — /ii). 
Consistency with the requirement of a crystalline con- 
figuration implies Tc ^ r < r*. Table 1 shows numerical 
values for a few favorable A^, where this process could 
be visible. Figure [2] shows our results for A^ = 12. Fully 
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FIG. 4: Finite temperature transition from the non-classical 
to the classical crystal for r = 10 and = 12. (a) and (b) 
are Monte Carlo snapshots for T — 0.00125eo and 0.005eo, 
respectively. 
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consistent with our estimates, we find that in the range 
3 < r < 15 the ground-state is the ring structure (4,4,4), 
[panel (b)], while the classical crystal configuration 
(3,3,6) sets in for r > 15, [panel (d)]. These two 
configurations have rotational symmetry C4 and C3, 
respectively. Actually, quantum mechanical fluctuations 
modify the ring structure (4,4,4) into two concentric 
rings with 4 and 8 particles, respectively. The radial 
density profiles for some r are shown in panel (d). Most 
interestingly, we find that in the intermediate range of 
T where the non-classical crystal is the ground-state, 
the classical ground-state can survive as a low-energy 
metastable configuration. It is then possible to induce 
a non-classical/classical-crystal transition by raising the 
temperature for fixed r. This is shown for = 12 and 
T = 10 in Fig. [31 Further melting into a featureless 
normal fluid occurs at even higher temperatures, which 
we find consistent with the classical melting ternperature 
for the homogeneous crystal configuration [13[ (not 
shown). Again, consistent with our semi-classical model, 
we found that the non-classical configuration (1,4) is not 
realized for N ~ 5, since ~ 1 < Tc. 

The mesoscopic phases above can be realized with po- 
lar molecules of current experimental interest. For exam- 
ple, for a moderate in-plane confinement uj/2tt = IkHz, 
ro = 0.4,0.8 and 1 urn for RbCs (d = 1.25Debye), LiCs 
{d = 5.5Debye), and SrO molecules {d = 8.9Debye), re- 
spectively, and r ~ 4,9 and 11. In addition, we pro- 
pose to observe the crystals using cold alkali atoms, by 
weakly dressing the groundstate of each atom with an 
excited Rydberg state with large dipole moment dr, in- 
duced by a strong DC electric field £bg directed per- 
pendicular to the 2D plane. For large interparticle sep- 
arations r ^ 100 nm, outside of the dipole-blockaded 
region [^, [l^ , the effective interatomic interaction is well 
described by the dipole-dipole potential with an effec- 
tive dipole-moment d ~ {n/A)'^dr, where fl and A are 
the laser Rabi frequency and detuning from resonance, 
respectively. The effective single-particle spontaneous 
emission rate is F ~ {U/A)'^Tr, with F^ the bare value. 
If one translates the latter into an effective heating rate, 
we expect the crystalline phases (r » 1) to be observable 
up to an average "melting" time T ~ rM/(Feft-Ei?), with 
Eji the photon recoil energy (in the tens of kHz), and 
Tm ~ O.leo the estimated melting temperature in the 
trap [l3|. For example, the groundstate of ^^Rb atoms 
can be weakly coupled to the Rydberg-statc {n = 20, 
ni = 19, 71.2 — 0) with bare dipole moment dr ~ 1.45kD 
and Tr/2Tr ^ lOOkHz. Here, ni and 712 are parabolic 
quantum numbers For a DC field Sue ~ 25kV/m 
and a coupling laser with n/2TT = 80MHz and A/2tt ^ 
IGHz, we obtain F/27r « 640Hz, Er « /i25kHz and 
d Ri 9.1Debye. For in-plane confinement tj/27r = IkHz, 
this gives ro w 1.07/im, r w 10 and T w 62.5/iS. To 
check the estimate for T we performed molecular dynam- 



ics simulations for the 2D system in the classical limit 
(t — > 00), where the spontaneous emission was simulated 
by applying random in-plane kicks to each particle with 
momentum (2mi?fl)^/^ at an average rate F. We found 
the average energy {E{t)) to increase linearly with time 
t. For the parameters above and A^ = 13 particles, the 
value {E(t) - E{0)) /2N sa 0.1 x eo is reached at t « 20/iS 
(roughly corresponding to a single kick) consistent with 
the above estimate for T. However, the actual lifetime 
(corresponding to the thermal equilibration of the sys- 
tem) is found numerically to be of the order of T ~ 200/iS, 
which is (much) larger than T [l^. 

Since the in-situ interparticle distances can be of the 
order of a /im or more, it may be possible to directly 
address single particles in situ, and thus image the spa- 
tial structure of the crystalline phases above using, e.g., 
tightly focused beams. Alternatively, we propose the fol- 
lowing method, which amounts to a cold-atom version of 
a magnifying lens. At a given time to the (DC or AC) 
fields inducing the dipole-dipole interactions are switched 
off, and the in-plane harmonic confinement is inverted in 
sign. Because of this inverted potential, each particle ex- 
periences a radial acceleration which depends on its spa- 
tial position at time to- After a certain time-of- flight, the 
particles can be, e.g, ionized and their positions recorded 
on a ion plate with unit efficiency, providing a magnified 
picture of the in-situ spatial configuration. 
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